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Abstract 

The method of dimensional recurrences proposed by one of the authors [l] [2] is apphed to the 
evaluation of the pentagon-type scalar integral with on-shell external legs and massless internal 
lines. For the first time, an analytic result valid for arbitrary space-time dimension d and five 
arbitrary kinematic variables is presented. An explicit expression in terms of the Appell hyperge- 
ometric function and the Gauss hypergeometric function 2F1, both admitting one-fold integral 
representations, is given. In the case when one kinematic variable vanishes, the integral reduces 
to a combination of Gauss hypergeometric functions 2Fi- For the case when one scalar invariant 
is large compared to the others, the asymptotic values of the integral in terms of Gauss hypergeo- 
metric functions 2^1 are presented in c? = 2 — 2e, 4 — 2e, and 6 — 2e dimensions. For multi-Regge 
kinematics, the asymptotic value of the integral in d = 4 — 2e dimensions is given in terms of the 
Appell function F3 and the Gauss hypergeometric hmction 2-^1. 
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1 Introduction 



Theoretical predictions for ongoing and future experiments at the CERN Large Hadron Colhder (LHC) 
and an International e'^e~ Linear Collider (ILC) must include high-precision radiative corrections. The 
complexity of the evaluation of such radiative corrections is related, in particular, to the difficulties 
in calculating integrals corresponding to Feynman diagrams with many external legs depending on 
many kinematic variables. Purely numerical evaluation of such integrals cannot provide sufficiently 
high precision within reasonable computer time. The evaluation of one-loop integrals corresponding 
to diagrams with two, three, and four external legs was studied in numerous publications. 

As for integrals associated with diagrams with five and more external legs, the situation is quite 
different. Not so many results for such integrals are available in the literature. Various authors 
[3] discussed the reduction of pentagon integrals to box integrals in space-time dimension d = 4. 
They were able to express these integrals as linear combinations of five different loop integrals with 
four external legs. Infrared divergences if any were supposed to be regulated by introducing a small 
fictitious mass. The representation of the dimensionally regularized pentagon integral in terms of 
box-type integrals was considered in Refs. [HIS]. However, in the calculation of multi-loop radiative 
corrections, higher orders in e = (4 — d)/2 are needed, and, therefore, one should extend such an 
expansion beyond the "box approximation." The first step in this direction was recently taken in 
Ref. [6], where an analytic result for the one- loop massless pentagon integral with on-shell external 
legs as well as several terms of its e expansion were presented in the limit of multi-Regge kinematics. 
However, no analytic results for arbitrary kinematics are available until now. Practically nothing 
is known about the analytic structure of on-shell pentagon integrals with unconstrained kinematics 
in arbitrary space-time dimension d. For very simplified kinematics, a pentagon-type integral for 
arbitrary value of d was given in Ref. [7] in terms of Euler gamma functions. A representation of the 
pentagon integral in terms of a four- fold Mellin-Barnes integral may be found in Ref. [8]. 

Important applications that require the evaluation of Feynman integrals with massless propagators 
include the study of jet production in QCD [9], which allows for a high-precision extraction of the 
strong-coupling constant a^, the investigation of the iterative structure of A/" = 4 supersymmetric 
Yang-Mills (SYM) amplitudes [10], and tests of the scattering-amplitude/Wilson-loop duality 

In this paper, we perform a first analytic study of the on-shell pentagon integral with massless 
internal lines and arbitrary kinematic invariants. We use the method of dimensional recurrences 
proposed in Refs. [BIS]) which was already applied to the calculation of one- and two- loop integrals in 
Refs. [2l[12l[13] and, quite recently, also to the calculation of three- and four-loop integrals in Ref. pH]. 

This paper is organized as follows. In Section 2, we introduce our notations, explain the recurrence 
relation for the pentagon-type integral with respect to the space-time dimension d, and we present 
a detailed derivation of its solution in Section 3. In Section 4, we consider a particular case of the 
on-shell pentagon integral with one vanishing kinematic variable and present an analytic expression 
in terms of the Gauss hypergeometric function 2^1- In Section 5, we specify the asymptotic values of 
the pentagon integral in d = 2 — 2e, d = 4 — 2e, and d = 6 — 2e space-time dimensions when one of the 
scalar invariants is much larger than the others. In Section 6, we present an analytic expression for the 
pentagon integral in the limit of multi-Regge kinematics in terms of the Appell function F3 and the 
Gauss hypergeomteric function 2F1. In the Conclusions, we summarize the accomplishments of the 
present paper and offer some perspectives for the application of the method of dimensional recurrences 
to six-point integrals. In Appendix A, we collect useful formulae for hypergeometric functions used in 
this paper. Appendix B contains intermediate results from Section 3. 
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2 Definitions and dimensional recurrences 

We consider the following integral with five massless propagators: 



4'^^('Sl2,S23,'S34,S45,Sl5; Sl3, ^U, ^24, ^25, S35) = / T^^IJ^: (2-1) 

i=i 

where 

D,={q-p,f +ie. (2.2) 

The pentagon diagram associated with the integral I^\{skr}) is depicted in Fig. 1. The labeling of 
the momenta in Fig. 1. corresponds to Eq. ()2.ip . The Lorentz invariants are defined as: 

Sij = Pij, Pij =Pi- Pj- (2.3) 

In the present paper, we take the squares of the external momenta to be vanishing, 

Sl2 = S23 = S34 = S45 = Sl5 = 0, (2.4) 

and work in the Euclidean region, Sjj < 0. In what follows, we only keep non-vanishing variables as 
arguments of ijf'^ and use the notation 

4'^Hsi3, Sii, S24, S25, S35) = /s"*^ (0, 0, 0, 0, 0; Sl3,Sl4,S24,S25,'S35), (2-5) 

for the on-shell case of Eq. (12. 4p . Due to the symmetry with respect to permutations of all propagators. 




the integral considered as a function of the kinematic variables Sij must fulfill the following relations: 

4'^^(S13, Sl4, S24, S25, = if^ (si3, S35, S25, ^24, S14) 

= /5'^^(S14,S13,S35,S25,S24) = ^''^ (^14, ^24, ^25 , ^35 , S13) 
= 4'^^('S24,S25,'535,Si3,Si4) = /s'^^ (S24, ^14, Sl3, 535, S25) 
= 4'^^(S25, 524,514, 513,535) = /g'^^ (525, 535, S13, S14, S24) 

= 4^^^(535,513,514,524,525) = 4*^^ (535, 525, 524, 5l4, 513). (2.6) 

2 



As was shown in Refs. [HIS], the integrals /^"^^^^ and /g*^^ fulfill the following relation: 

{d - 4)si3S24S35Sl4S25/l54'^^^^ (■^IS; ^14, S24, ^25, Sss) 

= 2S13S24S35S14S254'^^(S13, S14, §24, ^25, S35) + P^'^K (2-7) 



where 



k=l 

^\ = -S24S35('S24S35 " S^bSlA - S13S24 + S13S25 + S14S25), 
><^2 = -S35S14(S35S14 + 513^24 " ^24535 " S14S25 + S13S25), 
><3 = -S14S25(S14S25 " S35S14 " 513^25 + S24S35 + S13S24), 
>'^4 = S13S25(S13S24 " ^13525 " S24S35 " 535^14 + S14S25), 

>!^5 = S13S24(S13S25 " S13S24 " S35S14 + S24S35 " S14S25), (2-9) 

Jl = 4''^ (0,0,0,525; 524,535), J2 = 4"*^ (0,0,0, S13; Si4,S35), 
•^3 = 4^^^(0,0,0,524; S25,5l4), J4 = 4"*^ (0, 0, 0, S35 ; Sl3,S25), 

^5 = 4''^ (0,0,0,514; S24,S13), (2-10) 

r22 2 22222222 

/15 = [-524^35 ~ 2S24S13525 + •Sl3'S24 + •Sl3'S25 + •5l4'535 + 'Sl4'525 

Sl3S24S35'Sl4S25 

+2S13S24S35S14 - 2s24Sl4si5 - 2S13S35S24 + 2s24Sl3'S35'S25 + 2S24S35'Sl4S25 
-2S35S14S25 + 2S13S35S14S25 - 2S13S14S25 + 2S24S13'S14'S25] • (2-11) 

Here, I^'^ are integrals corresponding to Feynman diagrams with four external legs defined as 

M), . x_ [ d'^q 1 . ^ 

-'4 {Snj ■, Sjii, Ski, Sfih SjliSnk) — J £) D, ]J ,]J ' {■^■^■^J 

where Dj are defined in Eq. (j2.2p . An analytic expression for this integral with arbitrary kinematics 
in arbitrary space-time dimension d was recently obtained in Ref. |15] . 

3 Solution of the dimensional recurrence relation 

Redefining the integral 4^^^ as 



Ad)f N ^5 2 -(d) , . /Q 1 Q\ 

^5 (•S13, Sl4, S24, S25, S35j = , h (■Sl3, Sl4, ^24, ^25, S35), (3.13) 



we obtain the relation 

^5 (•Sl3,Sl4,S24,S25,S35j = ^5 (^13, Sl4, S24, ^25, S35j + Tj 

2S13S24S35S14S25 

which has the following solution: 



4n^l3,Sl4,S24,S25,S35)=n„((i) + ^- ^ ^ ' hi ' ' pi'^^^r^^) , (3.15) 



„ 2S13S24S35S14S25 
r=\j 
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where Ilm{d) is an arbitrary periodic function depending on the scalar invariants Sij and satisfying 
the condition 

U^{d + 2) =U^{d). (3.16) 

Another solution of Eq. (j3.14p reads: 

/r(si3, SU, S24, S25, ^35) = ^pid) - V ^ ^ ' hf' P^''+^^\ (3.17) 

~^ 2S13S24S35S14S25 

where 

Iip{d + 2) =Iip{d). (3.18) 

The correctness of both solutions, Eqs. (j3.15p and (|3.17p . may be easily verified by direct substitution 
into Eq. (|3.14p . Thus, for example, using Eqs. (j3.17p and (|3.18p . we have 

^ 2S13S24S35S14S25 

'd+2r-4'i 



00 p ( d+2r-A \ ^ 

7^1 2S13S24S35S14S25 



= 4''\^13,^14,S24,S25,S35) + ^ ' ^ 4 P^''\ (3.19) 

2S13S24S35S14S25 

in agreement with Eq. (I3.14p . Furthermore, the solution in the form of Eq. (I3.17P may be easily 
obtained from Eq. (j3.15p by adding to and subtracting from the expression on the right-hand side of 
Eq. (13.15P the sum 

y ^ ^ ^ h-f P^''-^'-^^- (3.20) 



_^ 2S13S24S35SI4S25 



In fact, the combination 



( d-2r-6 \ d ^ -00 p ( d-2r-6 > ^ 
^ \ 2 ) J-^^^ D(d-2r-2) , V- i ' -~ ' 



V 2 J 1^2-''^^ p(d-2r-2) _^ ^12; f^2-'"-'-p(d-2r-2) 

2si3S24S35'Sl4S25 ^ ^^^^ 2S13S24S35S14S25 ^ 



00 p / d-2r-6 \ d 

y ^ ) ^-r-lpid-2r-2) .321) 



2S13S24S35SI4S25 

is invariant with respect to the change d — t- dziz2l, where / is integer, so that this sum may be absorbed 
into the periodic constant that we denoted by Ilp{d). Changing the summation index in the remaining 
sum as r — >• —r — 1, we obtain Eq. (j3.17p . To obtain the solution of the difference equation (j3.14p in 
terms of convergent series, one may choose either Eq. (j3.15p or Eq. (j3.17p depending on the kinematics. 

The dependence of the arbitrary periodic functions Ilmid) and Ilp{d) on the scalar invariants Sij 
may be constructed from a system of differential equations which follows from the one for the integral 
I^^^ . For the integral Jg*^^ , we derive a system consisting of 5 differential equations of the form 

54^*^(513,514,524,525,535) r;\ 7 i?W 

S13S14S24S25S35 Sij ^ —— = [d-5) JkR' 



+ 



fc=l 

d dh^\ {d), . ,^ „„s 

2'5l35l4S24S25'S35'Sij^^ j ^ (Sl3 , ^M, ^24 , •S25 , •S35) , (3.22) 



where (i,j) = (1, 3), (1, 4), (2, 4), (2, 5), (3, 5) is not summed over, are the box integrals defined 
Hj 



in Eq. (I2.10p . and i?-^'' and are polynomials in Skr- To derive this system of equations, we use 
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the method proposed in Ref. [16]. Some details of the derivation are presented in Appendix B. The 
derivation of such equations is done with the help of the computer program package Maple. Explicit 
expressions for R^j'^ and are given in Eqs. (jl0.76p and (|10.77p in Appendix B, respectively. Using 

Eq. (j3.13p , one may obtain from Eq. (j3.22p a system of equations for the integral I^^^ . Substituting 
Eq. (IS.lSp into this system, we obtain the following system of equations for the periodic function 
Ilm{d) after a rather tedious calculation: 



CS13 '15'513S14S24S25S35 OSu ft5Sl3Sl4S24S25S35 



dUm{d) (j)24 , dUmid) 



S2A—^ = T ilm(aj, ^25^ = T o ^^rn\d), 

OS2A n.5'Sl3Sl4S24S25S35 C'S25 HsnSliS2iS2^S2,5 

sj-^ = limid). (3.23) 

CS35 n'5'Sl3Sl4S24S25S35 



These differential equations do not depend explicitly on d and are, thus, much simpler than those 
for the integral /g*^^ itself. The solution of this system of differential equation with respect to Sjj for 
Hrn{d) obtained with the help of the computer program package Maple reads: 



n„(d) = ,3 , (3.24) 

(Sl3'Sl4S35S24S25)'' 

where K{d) is an arbitrary periodic constant, 

K{d + 2) = K{d), (3.25) 

which is independent of the scalar invariants Sij. The system of differential equations for Tip{d) looks 
similar to Eq. (j3.23p . An arbitrary periodic constant K{d) may be determined from l'^^ calculated for 
some particular kinematics. Usually, setting some of the scalar invariants Sij to zero greatly simplifies 
the computation of the integral. But in our case, as may be seen from Eq. (j3.24p . K,{d) cannot be 
determined from l'^'^ calculated for such kinematics because the term proportional to K{d) drops out. 
For the same reason, one cannot use I^"^ calculated for kinematics with /15 = 0. Instead, we determine 
the periodic function Hrn{d) by comparing the limiting value of our analytic result for — )• 00 
with the analogous value obtained from the integral representation of Eq. (|9.63p by exploiting the 
steepest-descent method described in details in Ref. [12j. 

Without loss of generality, we henceforth assume the following hierarchy between the scalar invari- 
ants: 

- si3 > -si4, -S35 > -S25 > -•S24- (3.26) 
For the case when the scalar invariants Sij satisfy the conditions of Eq. (13.26P and additionally 

- si4 - S35 + si3 > 0, (3.27) 

the integrals I^f^ in Eq. (|2.10p may be written as 

0, 0, S25; S24, S35) = X^'^nS25, S25, ^SS, ^24) " ^^'^''(S24, S25, ^SS, •§24) " X^'^Hs35, ^25, ^35, S24), 
0, 0, S13; Si4, S35) = X^'^\si'i, Si3, Si4, S35) - X^'^Hs35, Sl3, S14, S35) - X^'^\sii, S13, S14, S35) 



2Tl{d - 3) ^[d) f -535514 



S35Sl4Sin^ ^ V'SSS + 514 - Sl3 

4'^^(0, 0, 0, S24; S25, S14) = ??^'^H^24, S24, S25, S14) - X^'^Ksu, S24, ^25, ^m) " X^'^ {S2^-, ^24, ^25, S14) 

_ 47r(d - 3) ^(d) / ■§25^14 

Sl4S25tan^ ^ V^25 + Sl4-S24 
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4'^''(0, 0, 0, S35; Si3, S25) = X^'^\S35, S35, Si3, S25) - V^'^KS25,S35, ^IS, S25) " X^'^HsiS, ^35, Sl3, S25) 



Ad) 



(d), 



Ad) I 



4'^^(0, 0, 0, S14; S24, S13) = X^'*''(S14, Sl4, S24, S13) " ^''"'^513, Sl4, S24, ^is) 

- ??^'^H'S24,Sl4,S24,Sl3), 



(3.28) 



where 



r/('^^(s,Si4,S24,Sl3) 



4(d-3)jf^(. 

(d - 4)si3S24 



■ 2-t^l 



2 
d _ 
2 



2; S(si3 + S24 - S14) 



1 



Sl3'S24 



X^'^^(s,Si4,S24,Sl3) 



-4(d-3)/f (. 



{d - 6)s(si3 + S24 - S14) 



1,3 



S13S24 



A-d. 
^ 2 ' 



S(S13 + S24 - S14) 



(3.29) 



and l2^\p'^) is the one-loop massless propagator-type mtegral 



ITT" 



d'^q 



- I 2^-- 

-vra {—p )2 



dj'i J q^[q — 



2^-3r(^)sin , 



d • 



(3.30) 



An expHcit derivation of these results using the method of dimensional recurrences may be found in 
Ref. [15]. Results for these integrals were also obtained in Ref. [1] using a different method. 
To evaluate I^^'^ , we use the solution in the form of Eq. ()3.15p . This solution may be used if 



< 1. 



(3.31) 



As we shall see later, the quantities ^/{h^Sij) emerge as expansion parameters in the resulting hyperge- 
ometric series. If condition (j3.3ip is fulfilled, then we may obtain a convergent series using Eq. p.l5p . 
If this condition is not fulfilled, then we may use Eq. (I3.17p . In this case, the inverse quantities, 
h5Sij/4:, will be the expansion parameters in the resulting hypergeometric series. 

In the following, we obtain an analytic result assuming that Eq. (j3.3ip is fulfilled for all scalar 
invariants Sij. If this condition is not fulfilled for a particular term in P^^\ then an analytic result 
may be obtained by performing analytic continuations of the hypergeometric functions in the final 
result. Another possibility to obtain the analytic result in this case is to repeat the calculation using 
the solution of the form of Eq. ()3.17p . It should be noted that, if Eq. (|3.27p is not satisfied, then the 
arguments of all hypergeometric functions generated by the integral ^4'^^ (0, 0, 0, S13; su, S35) exceed 
unity. Analytic continuation of the result for this integral given in Eq. (j3.28p yields 



4'^^(0, 0, 0, S13; Si4, S35) = r/^'*^(si3, Sl3, Si4, S35) - r?^'*''(s35, Si3, Si4, S35) - ??^''^(S14, S13, S14, S35). (3.32) 

Adopting Eq. ()3.3ip . exploiting Eq. ()3.28p and, if — S14 — S35 -|- S13 < 0, also Eq. (j3.32p for J2, we 
obtain a result in which each term is real. 

The infinite sums resulting from Eq. (j3.15p may be written in terms of known hypergeometric 
functions. As is evident from explicit expressions for the I^'^^ integrals, we must compute three 
different types of sums. The first one is related to the function x^'^K 



Ad) 



id) I 



Ad) I 



r=0 



-4hi 



d-2r-6 
2 

-1 00 



h 2 
"-5 



s Q 

r=0 

A{d-3){d-5)r 



-4) 



3) {d-2r- 5) 



{d - 2r 



2F1 



Q s^ix-1) 



4'\s) F3 



1,1, 



l,r + 4 
r + 5 - 

10 -d 



_ d. ■ 

2 ' ^ 

d . ^ 

2 ' 



2 ' s/i5 ' a; — 1 



(3.33) 
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where 



Q = c + b — a, 



(3.34) 



The result of the summation of a term including the r]^'^^ function reads: 



r=0 



d-2r-6 



h 2 

"-5 



-r-1 



8(d-3)(d-5) ^f-i.^ /^d - 6 
6cs ^ 



8^) Vs/^sy {d-2r-Q) 



l,^-r; s{c + b-a) 

2 



r : 
r ; 



be 



(3.35) 



The simplest type of infinite series originates from the two terms in Eq. (|3.28p without 2^1 functions. 
The contribution arising from the term proportional to X2 reads: 



d - 2r - 6 



\ J-r-l (-2^)(d-2r-5) ^(d^2r-2) 
) % . ^(<i-2r-2) ■'2 

/ s,s5Si4 sm^ o i 



r=0 V - / 535514 

^ 87r(d - 3)(d - 5) (514 + ^35 - Si3) ^(d) 
(d-6)sin^' 



-S35S14 



S35 + Sm - Sl3 



„2 2 
■'35*14 

7-d . 



-S35'Sl4 



S35 + Sl3 - Sl3 



1, V; 4(g35 + 514 - gl3) 
"2" ' -535514^5 



(3.36) 



and the one related to X3 reads: 



Er 

r=0 



d-2r-6\ ,f-r-i (-47r)(d-2r-5) 



"5 



~2r-2) 



, 7r(d-2r-2) 2 

S14S25 tan 2 ^ 



S14S25 



S25 + Sl4 - S24 



167r(d-3)(d-5) (si4 + ^25 - ^24) ^i-ip 
(d-6)tan -^^ "2 „2 5 



X 2i^l 



2 *35'5i4 
^ • 4(S14 + S25 - S24) 



(d) 



S14S25 



S25 + Su - S24 



(3.37) 



—2~ ; S25Sl4^5 

Inserting the sums of Eqs. p.33p . (I3.35p . (j3.36p . and (I3.37P in Eq. (I3.15p . we obtain the following 
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expression for the integral /g*^^ : 

4'^^(S13,S14,S24,S25,-S35) - 



5-d 



Sin^ ^ \/-«13Sl4S35S24S25 



+ 



^2(535 + gi4 - ■813) (c? - 8)TrK ^(d) 
2h,sl,sl, sinf 



-S35S14 



S35 + Sl4 - Sis 



2F1 



1,^; 4(535 + gi4 - gi3) 

; 535*14^5 



su + S25 - S24 



+ ^^^^4'^)(.24) f (^35 + ^24-. 25) \ 

^5*24 [535 S24 V^5S24 S35 / 



1, Y ' 

—2~ '■> S25Si4/i5 



•''■21) 



X3 



S14S25 



/ 4 g24(gl4 + S25 - S24:) \ _|_ 
1^5^24' 



S14S25 



S13S24 



/ 4 Sl3 + S24 - Su 
\h5S24 



Sl3 



K 



/i5Sl3 



4'^ (^13) 



X2 



S35S14 



+ 



X4 



(Sl3 - •S35)(S13 - su) \h^Slz' (si3 " S35)(si3 " Su) , 
$W ( -"25 ^ ^ ^ ^5 ^ ^(,) 4 -524 



(Sl3 - S35)S13 V^SSls' Sl3 - S35/ (si3 " Sl4)si3 \ ^5Sl3 ' ^13 " S14 / . 



^^^^^ 



X5 



X2 



S35 



(Si3 - S14)S14 V ^5^14 ' Sl3 - Sl4 / (si4 - S24)si4 



X3 



-S25 



/15S14' Si4 — S24 



— $W I ^ 

(Si4 - S24)(S13 - S14) \h^Su ' (Sl4 - S24)(si3 " Su) 



( 4 



S13S24 



K 



h^S2h 



4'\S2,) 



Xl 



( 



-S24S35 



(S25 - S35)(S24 - S25) \h5S25 ' (^25 " S35)(s24 - S25) 



^ Md - 8) 
S25S13 



(^5525' 



S25 + Si3 - S35 
Sl3 



+ 



X3 



(S24 - S25)S25 



Sl4 



h5S35 



Xl 



S24 



(S25 - 535)535 V^5535 S25--535/ (^13 " S35)s35 



+ 



X2 



^5^25 ' S24 — S25 



Si4 



X4 



.$(<i) 



S25S13 



where 



(S25 - •S35)(S13 - S35) V^5S35' (^25 - 535) (si3 - •S35)/. 

-8{d-3){d-5) 



K 



{d -6){d- 8)si3Si4S24S25S35 ' 

00 (z-,) 



\^5S35 S13-S35. 

(3.38) 
(3.39) 



J / 10-d \ 
V 2 A 



r + 5 - f ; ^ 



1,1, 



7 - d 10 - d 



,1, 



(3.40) 



00 fj-d'j 



r=0 



W 



(8^)^ (d-2r-6) 



2-t'l 



00 



1 £b6 
-•-50 ' , 



2 

(i-4 



— r : 



1 ^ / d— 6 \ 



1 6^ _ 7 . 

' 2 ' 2 ' 
8-d 8-d _ J . 
2 ' 2 ^ ' 



(3.41) 



and the quantities Xj are defined in Eq. ([IJ]). The functions and ^^'^^ are related as 



{l-z){d- 



1-z 



7r{-zf 



2sin^ 



2^1 



1 LzA ■ 

8-d . 



2 ' 



(3.42) 



The Appell function admits the following one-fold integral representation: 



1,1, 



7 



d 10 



■.w,z 



r(i)r(^) 



:i 



5-d 

n) 2 



arcsm vnz 



1 — w + uiy VI — 



du. 



(3.43) 



It is interesting to note that changing the space-time dimension d in Eq. (|3.43|) by one unit, 
d ^ d + 1, we obtain the Appell function i<3, which we already encountered in the calculation of the 
one-loop master integral entering the calculation of radiative corrections to Bhabha scattering [17]. 
As a matter of fact, we observed in Ref. [ITJ that, at the one- loop level, the set of hypergeometric 
functions appearing in the results for n-point integrals with massive propagators also appear (up 
to the change d — )• d + \) in the calculation of (n -|- l)-point integrals with all propagators being 
massless. For example, the one-loop propagator integral with two different masses is expressible in 
terms of two Gauss hypergeometric functions 2-^^1(1, {d — l)/2; d/2] z) with different arguments, while 
the result for the one-loop vertex integral with massless propagators is expressible in terms of the 
2-Fi(l, {d — 2)/2; {d — l)/2;z) function with different arguments. The result for the one-loop vertex 
integral with arbitrary masses and external momenta [T2] is expressible in terms of the Appell function 
Fi((d - 2)/2, l,l/2,d/2;x,2/) and the Gauss hypergeometric function 2i^i(l, (c? - 2)/2; (d - l)/2;z) 
with different arguments, while the result for the box integral with massless propagators and arbitrary 
external momenta is expressible in terms of the same functions up to the shift d — )• d+l, as was observed 
in Ref. [15j . We stress that the number of terms with hypergeometric functions appearing in the 
expressions for one-loop integrals with massive propagators and the arguments of these hypergeometric 
functions are different from those for their counterparts with massless propagators, but the sets of 
hypergeometric functions are the same up to the shift in d mentioned above. 

The analytic continuation of the result in Eq. (j3.38p amounts to the analytic continuation of the 
hypergeometric functions -F3 and 2-^1 and the factors in front of these functions. From Eq. (I3.38p . we 
may obtain the value of the integral in any region by using the usual ie prescription and observing that 

is manifestly real in the Euclidean region where all scalar invariants sij are negative. The analytic 
continuation of the Gauss hypergeometric function 2F1 is well understood |18] . Several useful formulae 
are given in Appendix A. The analytic continuation of the Appell function may be obtained from 
the series representation of Eq. (|3.40p by the analytic continuation of the 2F1 function. We notice that 
the well-known formula for the analytic continuation of the F^ function [18] in terms of the Appell 
function F2 is not applicable if both arguments of the F^ function in Eq. ()3.40p are large. In this case, 
one may proceed by analytically continuing the functions 2F1 and 3F2 in Eqs. (j3.40p and (I3.4ip . For 
example, the analytic continuation of the 3-F2 function in Eq. (j3.4ip is 



3-f^2 



1, 



7-d 6-d 
2 ' 2 



^ 2' ^ 2 



I; 
-I: 



W 



(6 - d){d-Q + 2l) 
w{d-5){d-4: + 2l) 



3F2 



1 £b4 dzd + 7 ■ 1 
' 2 ' 2 ^ ' 
d-3 d-2 I ; . ~ 



r(«^)(d + 2/-6)r(^) ^ 



3 . - 

2 ' W 



1 + 



2 sin^ r (^) l\ ^ ' ^ ' 

Substituting this formula into Eq. ()3.4ip . we obtain the following representation for the H^'^\w,z) 
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function for \w\ > 1: 



(6-d) 



-1) 



w 



2 sin^ r (^) VT 



wz 



\J vr {wz — 1) 



r 4 



d\^/d-5\ l{l-wz) 

- r arcsm \ / — ; , 

2) \ 2 ) Vwil-z) 



(3.45) 



where 



1;2;1 




2 ■ 2 

d-1 . d-2 
2 ■ 2 



, 1; 1; 



oo oo ( d-3 \ ( d-3 \ 
Sr^Sr^ \ 2 )r+s \ 2 )r r s 

r=0 s=0 \ 2 Jr+s \ 2 Jr 



CO (d^\ (d-3\ 
V ^ 2 )r r p 

/ . (d-2\ (d~l\ ^ 1 

r=0 V 2 )r\ 2 )r 



d-3 



2 

d-1 



r ; 



+ r : 



(3.46) 



and F^'i-Q is the Kampe de Feriet function [TU]. To analytically continue the H^'^\w,z) function into 
the region w,z > 1, one may use Eq. (j3.45p . The analytic continuation of the (p^'^') {w , z) function into 
the region |z| > 1 may be obtained by the analytic continuation of the 2-^1 function in Eq. (j3.46p using 
Eq. (I9.66P from Appendix A. 

There are several relations for the function which may be useful for its analytic continuation 
and performing its £ expansion. We present here two formulae for the analytic continuation of the -F3 
function with large first argument x and y < 0. One such relation follows from the results given in 
Ref. [17] and reads: 



-^3 ( l,l,-^Tr-'l'— TT— '^'^ 



+ 



2 

id -4) 
x{d-l)' 



{d-6){d-8) 



F3 7T,1>1, 



{d — 3){d — 5){xy — y — x)x 
d-2 d + l y 1 



2-ri 



y 



1 i- 

' 2 ' 
d-1 . ' , 

-5- , x + y - xy 



/vr 



d 



ri^-^]Vi-x{-x) 



2 X + y — xy X 

d-S 



2F1 



1,1; y{x-l) 

3 . 

2 ' 



X 



(d- 



(d-5) x(l-y) 



iFi 



d-4, 
2 ' 

2 ' 



x{y-l) 



iFi 



1,3 



2 ' 

d . 

2 ' 



y 



y-i 



(3.47) 



The -F3 function with appropriate parameters admits the one-fold integral representation 



^ , 1 ^ d-2 d+l 
F3 ( 2'1'1'^— ; ^'^ 



r(^; 



^ r (^)w Jq 1-z + vz 



(l-v) 



d-4 



'WV 



1 



dv, 



'WV 



(3.48) 



which may also be used for the analytic continuation. Another formula for the analytic continuation 
into the region of large x and y < connects the function and the Horn function H2, 



7-d 10 -d 
F3 1, 1, -T^, 1, — 7: — ;x,y 



d-5 



r 



i8-d) 
x{d-5) 

10 -d 



Ho 



d-6 d-3 1 



(-x) 2 ^/l- X 2F1 



1,1; y{x-l) 

3 . 

2 ' 



X 



(3.49) 



We checked Eq. ()3.38p by comparing its numerical values calculated for some specific values of d and 
Sij satisfying the condition of Eq. (|3.26p with the result of the direct numerical evaluation of the four- 
fold integral representation given by Eq. (|9.63p in Appendix A using the computer algebra program 
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package Maple. In all cases, we found perfect agreement to all valid digits of the numerical results 
(usually about 12). For example, we performed a numerical calculation for d = 21/2 and d = W — 2e, 
where e = 1/10000, setting S13 = -1000 + 6, su = -65 + 6, S35 = -190 + 6, S24 = -1/20 + 6, and 
= -1/10 + 5 with 6 = ilO-20. 



4 Case of one vanishing variable 



In practical applications, such as the analytic continuation of the pentagon integral with the help 
of functional equations |17| I20| . we need the value of the integral /g'^^ with one or more vanishing 



kinematic variables. As was shown in Ref. [20], the one- loop vertex integral with arbitrary masses and 
external momenta may be expressed in terms of a vertex integral, in which two masses and the square 
of one external momentum vanish, with the help of functional equations. The formulae for the analytic 
continuation of the vertex integral with arbitrary masses contain integrals with all masses vanishing. 
The relation between the master integrals of Bhabha scattering and heavy-quark production derived 
in Ref. [T7] includes an integral with massless propagators. In all these cases, relations connecting 
integrals with different kinematical variables include integrals with simpler kinematics, in which either 
some masses or some squared momenta are equal to zero. We observed a rather similar situation 
for the pentagon integral. Functional equations relevant for the analytic continuation of a considered 
integral include the value of this integral with one scalar invariant taken to be zero. For this reason, we 
wish to present here the value of the integral /g'^'* with one scalar invariant taken to be zero. Due to the 
symmetry (|2.6p of this integral, it is sufficient to consider the case when S24 = 0. A detailed discussion 
of functional equations for pentagon-type integrals will be presented in a separate publication [15] . 
For S24 = 0, we obtain the following result: 

4'^''(si3,si4, 0,825,535) 



-8(d-3)(d-5) 



+ 



{d - 6){d - 8)(si4S35 - S14S25 + S13S25) 
S35S14 



irjd - 8)(si4 + S25) j(d) f \ 

Sl4S25tan^ ^ V^14 + '525y 



(Sl3 - S35)(S13 - SiAjSi-i 



(^13)2^1 



1,1; 

d . 



2 

S35'Sl4 



5 - 2 ; (Sl3 - S35)(S13 - Sli) 



^25 j{d)f^ ^ p 

H -^2 lSl3j2-C'l 



SlSl^ia - S35j 



■S35 j{d) 



Sl4(,Sl3 - SlA, 



1, 1 ; •S25 



5 — 2 ! ■535 — S13 
1,1; •S35 

r d . 

^ ~ 2 ' ^13 — Si4 



— I2 (514)2-^1 



'14 



1, 1 ; S25 



2 ' 



Sl4 



- — -'2 {^25)2^1 
^25 



1,1; si4 



2 ' 



S25 



id -8) id). . J,' ' " 

-'2 \S25)2^] 



{d - 6)S25' 



l,f-3; 525 + si3 - S35 
f — 2 ; S13 



-'2 lS35j2-t'l 



+ 



S35[S13 - S35 

S13S25 



1,1; Si4 



5 — 2 ; ■Sis — 



535(525 - S35)(S13 - S35; 



1^(535)2^1 



1,1; 

d 



S13S25 



5-2; (525 - S35)(S13 - 535) 



+ 



^(■Sl3 - ^35 - Sl4)(^ - 8) ^(d) f S35 



2sin^ S14S35 ^ VS13 - S35 - 



(4.50) 



We observe that this result is significantly simplified compared to Eq. ()3.38p . In fact, only 2F1 functions 
remain. 
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5 Asymptotic values of the integral Jg*^^ for |si3| — )- oo 

If the scalar invariants fulfill the hierarchy of Eq. (j3.26p with |si3| being much larger than the other 
scalar invariants and S24 / S25, then we have 



[S24 - S25) 



-Sl3 



S24'S35Sl4S25 

2(si4S24S35 - SuS 



+ 



2 

25 



■§24^*35 + S24S35'S25 + S14S25S35 + S14S24.S25) 



•S24S35"Sl4'S25 



+ 0{ — 

Sl3 



(5.51) 



id) 



As one can see from the explicit expression of Eq. ()3.38p . the first argument of the two functions ^^'^'^ 
and H^'^'^ is always proportional to l/|si3|, so that it is sufficient to keep only the first terms of the 
series in Eqs. (|3.40p and (j3.4ip in order to find the asymptotically leading term of the integral 
In d = 4 — 2e dimensions, this value reads: 

H (Sl3,Sl4,S24,S25,S35) 

_ 7r(l - 2e) (S25 + S24) ,-{4~2e) / S35S14 



[si3|->-oo 



p 

smvre S14S35S24S25 ^ 

7r2r(l+e) (-Sl3r[(s24 



Sl3 



^25)^]^+^ 



sin^vre 



Ind 



(S14S35S24S25)^+^ 

2 — 2e dimensions, the result is 



+ 



{S25 - S24) 

S25{S25 + S24, 
1 



■ 2F1 



1, 1 - e; £24 
1+e; S25 



\S13 



l-e 



(5.52) 



r(2-2£)/ s 
^5 [Sl3,Su,S24: 825,835) 



^ 4TT{l-4e^){3 + 2e)is25 + 824)si3 

|si3|-s>oo 8l^sl^{s24: - 525)^(2 + e) sinvre 



X /: 



{4~2e) / S35S14 
813 



-513(^24 - 525 ) 
Sl4'S35'S24S25 



2F1 



1 



- +£■ 

' 2 ' 

3 + e; 



-4S24S25 
{824 - S25)' 

7r2r(2 + e) 



\/-'5l3Sl4S24S25S35 sin TT£ 



+ 



I •5131 



(5.53) 



We observe that, when one of the scalar invariants is large compared to the others, the hypergeometric 
functions of two variables collapse to hypergeometric function of just one variable. 

In d = 6 — 2e dimensions, the leading large-|si3| term of the integral li'^^ is somewhat more 
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complicated, being 
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S14S25 
(4-2£)/_.^ 7-(4-2£) 



+ 7^^^) (524) In 
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^14 r(4-2e) 
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1,1 +£; S24S35 
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S24 



(1 + e) (535 + 524 - S25) 



4'"'^^ (^35)2^1 



S25 7-(4-2£)/ N p 
--'2 (•S35j2-ri 



1,1 +£; 525 

2 + £; 535 
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where ^(x) is the logarithmic derivative of the Euler V function, 

V'(x) 



|S13 



|2-£ 



r(x)' 



and 7 is the Euler-Mascheroni constant. 



(5.54) 



(5.55) 



6 Limit of multi-Regge kinematics 

In this section, we consider the pentagon integral I^"* with the somewhat peculiar kinematics, 

- 513 > -514, -535 > -525, -524- (6.56) 

In addition, we assume that the following relation among the scalar invariants Sij holds: 

514S35 ~ 513524 ~ 513525- (6.57) 

Equivalently, one may define such an ordering by introducing the scaling parameter A, putting 

513 513, 514 A5i4, 535 A535, 525 ^ A^525, 524 ^ A^524, (6.58) 

and taking the limit A ^ 0. In fact, this arrangement corresponds to the multi-Regge kinematics. 
Recently, the pentagon integral in d = 6 — 2£ dimensions in the multi-Regge kinematics defined by 
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Eq. (j6.58p was considered in Ref . [B] . The integral /^'^^ in d = 6 — 2e dimensions is needed to determine 
the e expansion of the integral /g*^^ in d = 4 — 2e dimensions. 

Here, we consider the integral /g*^^ in the multi-Regge-kinematics limit of Eq. ()6.56p directly in 
d = 4 — 2e dimensions. To obtain the leading contribution in this limit, we perform the scaling of 
Eq. (I6.58P and retain the leading terms in the limit A — )• 0. The most divergent terms in Eq. ()3.38p are 
of the order In A/A^'*'^'^. As will be seen later, such terms cancel out in the sum, so that the leading 
asymptotic terms are of order l/A^^^"^. In order to find the leading terms in the limit A — >• 0, we must 
know the asymptotic behavior of the H^'^\w, z) function for z — >• 1 and that of the <^'^'^\w^ z) function 
for z — )• oo. The leading and subleading terms of these functions in the respective limits may be 
obtained from the series representations of Eqs. (|3.4ip and (j3.40p by retaining the first leading terms 
of the expansions of the 2F1 functions. With the help of Eq. (I9.66P in Appendix A, we obtain 
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+ 0((l-2)ln(l-z)). (6.59) 



0=0 



In a similar way, the asymptotic behavior of the ^^'^'^ function for 2: — )■ 00 may be derived. Keeping 
the logarithmic and constant terms of Eq. ()9.68p in Appendix A, we arrive at the following result: 



w,z, 
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+ O . (6.60) 



Using Eqs. ()6.59p and ()6.60p . and retaining only terms contributing at orders InA/A^"^^^ and l/A^"*"^^, 
we find from Eq. (|3.38p in the limit of multi-Regge kinematics the leading term. 



where 
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(6.61) 



2S24S?3S25 + 2S13S35S14S25 + 2S13S24S35S14 



Sl3^S24S35Sl4S25 



(6.62) 
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As mentioned above, the terms with In A cancel, as well as those involving the function tp{x), 7 and 
the derivative of the 3F2 function with respect to a parameter. 

The value of the integral ijf"^ for d = 6 — 2e in the limit of multi-Regge kinematics may be derived 
from Eq. (|2.7p using Eq. (|6.6ip and the asymptotic values of integrals given in Eq. (|3.28p . We just 
note that it is again expressed in terms of -F3 and 2F1 functions because the integrals add only 2-^1 
functions. An analytic expression for the integral I^'^^ in d = 6 — 2e dimensions was recently obtained 
in Ref . [6] . The result of Ref. [6] is given in terms of derivatives of the Kampe de Feriet function [19j 
with respect to a parameter. For a direct comparison of our result with that of Ref. [6], one needs an 
expression of the function in terms of derivatives of the Kampe de Feriet function, which, to our 
knowledge, is not currently available. 



7 Conclusions 

In this paper, we evaluated the one-loop scalar pentagon integral in arbitrary space-time dimension 
d with on-shell external legs, massless internal lines, and otherwise arbitrary scalar invariants. Ex- 
ploiting the method of dimensional recurrences, we obtained a result in terms of the hypergeometric 
functions -F3 and 2F1. In our case, both functions admit one- fold integral representations suitable for 
e expansions. Using the methods of Ref. [6j, the on-shell pentagon integral may be represented in 
terms of four- fold hypergeometric series, while the method of dimensional recurrences advocated here 
just yields two-fold series. 

The method of dimensional recurrences may also be applied to the evaluation of the hexagon 
integral, which is needed for the calculation of the 0{e) contributions to one-loop maximally-helicity- 
violating amplitudes at one loop in = 4 SYM theory. To simplify derivations in this case, one may 
start from the dimensional recurrences written in the limit of multi-Regge kinematics. In our opinion, 
the method of dimensional recurrences is quite efficient to go beyond the "box approximation" for the 
n-point one-loop integrals. 

The £ expansion of our results in Eqs. (j3.38p and (I6.6ip will be presented in a future publication 

m- 

We expect that the results presented in this paper and our forthcoming one [15] may be conveniently 
incorporated in program packages for the automated analytic computation of one-loop integrals in 
massless theories, similar to the package for numerical calculations presented in Ref. [21]. 
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9 Appendix A 

9.1 Integral representation of I^'^^ 

For numerical checks of the result for the integral I^^\sis, s^, S35) in Eq. ()3.38p . we use the 

Feynman parameterization 

4'^^(si3,si4,S24,S25,S35) = -r ^5 - j ... j dxidx2dxsdxA xfxlx^ ^, (9.63) 
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where 



= XlX2 [(1 - X3){1 - Xl)si3 + XlX3X4{l - X2)s25 + X3X4{1 - Xl)s35 
+ XiX3{l - X4){1 - X2)S24 + 0:1X2X3(1 - X4)(l - X3)si4] . 



(9.64) 
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The integral representation of the -F3 function used for the derivation of the one-fold integral repre- 
sentations of the $ function reads: 



F3(a,a',/3,/3',7;x,y) 



r(7) 



1^7-/3-1(1 



r(7 - /3)r(/3) 70 {i-x + uxY 

This integral representation follows from Eq. (20) in Ref. |22j . 



7-/3; 



uy 



du. 



(9.71) 



10 Appendix B 



In this appendix, we present the system of differential equations for the integral l'^^ to be used for 
the derivation of the differential equations for the periodic constants Pp{d) and Pm{d). To obtain this 
system of differential equations for we exploit the method described in Ref. [16) . According to this 
method, derivatives with respect to Sij may be expressed in terms of integrals with shifted space-time 
dimensions. Explicit expressions for such derivatives may be derived from the integral representation 
with a parameters. 
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where 



(10.72) 



Q = aiassis + 0:104514 + 0204824 + 020:5825 + 0305835, 
D = Ol + 02 + "3 + "4 + 05. 



From Eq. (|10.72l) . it follows that 
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where 



^d+2 p 



d+2 



D1D3 

d'^+^q P 
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^d+2q p 

i-K^ D2D4 ' 
d'^+^q P 

iTT^ D2D5 
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d+2 
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D3D5 



(10.73) 



(10.74) 



(10.75) 



D1D2D3D4D5' 

and Dj are defined in Eq. (j2.2p . In order to reduce the {d + 2)-dimensional integrals on the right- 
hand sides of the relations in Eq. pu.74|) to a set of basic integrals, we use recurrence relations |23j . 
All calculations are performed with the help of computer program package Maple. The resulting 5 
differential equations for the integral ijf"^ all have the form of Eq. (j3.22p . The polynomials and 
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(f)ij occuring therein are found to be 



^13 = (-^25 + 524 + S35)(S14S25 - SSS^U - S13S25 + S24S35 + ■S13S24), 

(2) / 2 2 2 

^13 = -(^14*25 - S14S35 + S35S14S25 - S13S14S25 + 2S13S35S14 - S14S35 + S13S24S14 

+ S24S35S14 + S25S13S35 + 524535 - S24S13535), 

(3) 

RlS = -(-S24S35 + S35S14 + SI3S24 - SI3S25 - S14S25)(S14 + S25 - S24), 

-^if = ■Si3(-Si4S25 + S35S14 - 2^25 + 2524*25 " S13S25 + 2525*35 + S13S24 " -524535), 

-^13^ = -Sl3 (■514*25 - S35S14 - 2S14S24 " S13S25 " 2524*25 + -513524 + S24S35 + 2s|4), 

-^14 = (.513.525 - .513.524 + .524*35 + ^35.514 - .514.525) (.524 + .535 " .525), 

Ru = *14(*14525 - 524*35 + 2S13S35 - 535*14 " *13*25 + *13*24 + 2^25*35 " 2*15), 

Ru = S14(S35*14 - *14*25 " 2*^5 + 2s24*25 " *13*25 + 2S25*35 + *13*24 - 524*35), 

Ru = (*35 - *13 - *25) (535*14 - 524535 + 513S24 - 513S25 - 514525), 

Ru = (513524 - 524513535 + 524535514 " 2si3S245l4 - 524513525 " 535514 
+ 513S24 - 524514525 + 5l3Sl4525 " 513S35S14 - *i3*25), 

-R24 = 524(514525 - 524535 + 2S13S35 - S35S14 - S13S25 + 513S24 + 2^25535 " 25I5), 
(2) 

R24 = (-513525 + 513524 - 524535 - 535514 + 5l4525)(-535 - 514 + 513), 

R^4 = -(-5^525 + 5^4^35 - S14S25 + 513S14S25 + 535514525 

+ 2S245l4525 + 5l3S245l4 - 524535514 + 5135^5 + 524535525 - 524513525), 

-R24 = (513 + 525 - 535)(S35S14 + 513*24 - 524*35 - 514525 + 513525), 

-^24^ = -524(-524535 + 513S24 - 2si3Si4 - 2S13S35 + S13S25 + 2*?3 - S14S25 + *35*14), 

-R25 = -(-*24*35 + 2524*35*25 - S24513S25 + 525*13*35 

+ *24*13*35 + *13*24 + *24*14S25 - *35*24 + *14*35 " *35*14*25 + *24*35*14), 

(2) 

R25 = (-*13*25 + *13*24 + *35*14 - *24*35 + *14*25)(*35 + *14 - *13), 

-^25^ = *25(-2*i4 - *14*25 + 2si3Si4 - *35*14 + 2si4S24 + *13*25 - *13*24 + *24*35), 

-^25^ = -*25(-*24*35 + *13*24 - 2si3Si4 - 2S13S35 + S13S25 + 2*13 - *14*25 + *35*14), 

-^25^ = (*13 + *24 - *14)(*35*14 + *13*24 - *24*35 - *14*25 + *13*25), 

R35 = -*35(*14*25 - *35*14 - 2si4S24 - *13*25 - 2s24*25 + *13*24 + *24*35 + 2*14), 

-^35^ = *35(-2*i4 - *14*25 + 2si3Si4 - S35S14 + 2si4S24 + *13*25 - *13*24 + *24*35), 

-^35 = (*14 + *25 - *24)(-*13*25 + *13*24 + *35*14 - *24*35 + *14*25), 

R^35 = -(*14*25 - *35*14*25 + *13*14*25 + *13*35*14 - *13*25 - *13*25 

+ *24*35*25 + *24*13*25 + 2s25*13*35 - *24*13*35 + *13*24), 

-^35^ = (*14*25 - *35*14 - *13*25 + *24*35 + *13*24)(*13 + *24 - *14), (10.76) 
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and 



013 = -35^48^5 + 6S25S14(S14 - ^24)535 - S14S25(S24 - S25)S13 

- 3(si4 - 524)^^35 + 2(S24 " 825)'^ S^^ + (^24 - S14S25 " ^14524 " S24S25)'Sl3S35, 
(plA = -3S13S25 + (S35 - ■S25S35 - S13S35 - S13'S25)S14S24 

+ S13S25(S25 - S35)'Sl4 + 6si3S25(si3 " S35)s24 " 3(si3 - S35)^S24 + 2(s25 " S35)^Sl4> 
(j)24 = -3Si4s|5 - 6S35S14(S13 - Sl4)s25 " S35'Sl4(si3 " 835)524 

+ {^13 - S35'Sl4 - Sl3'Sl4 " Sl3'S35)'S24S25 + 2(si3 " S35)^'S24 - 3(si3 - Sl4)^S25' 
025 = -3Si4s|5 - S35Si4(si3 - Si4)s25 - 6535514(813 - 835)524 

+ (•^IS - S35S14 - S13S14 - Sl3S35)'S24S25 " 3(si3 - 535)^524 + 2(^13 - Sl4)^S25; 
035 = -3SmS25 + ■S25S14(S14 " ^24)835 " 6si4S25(s24 " 825)^13 + 2(si4 - ^24)^535 

- 3(S24 - •525)^'5i3 + (S24 " S14S25 " •S14S24 " •S24S25)'Sl3'S35, (10.77) 

respectively. 
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